Draft version November 12, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



PERIOD-LUMINOSITY RELATIONS DERIVED FROM THE OGLE-III FUNDAMENTAL MODE CEPHEIDS 

Chow-Choong Ngeow 

Department of Astronomy, University of Illinois, Urbana-Champaign, IL 61801 



Shashi M. Kanbur 
Department of Physics, State University of New York at Oswego, Oswego, NY 13126 



HiLDING R. NeILSON 

Department of Astronomy & Astrophysics, University of Toronto, Toronto, ON, Canada M5S 3H4 



A. Nanthakumar 

Department of Mathematics, State University of New York at Oswego, Oswego, NY 13126 

AND 

John Buonaccorsi 

Department of Mathematics & Statistics, University of Massachusetts, Amherst, MA 01003 
Draft version November 12, 2008 



ABSTRACT 

In this Paper, we have derived Cepheid period-luminosity (P-L) relations for the Large Magellanic 
Cloud (LMC) fundamental mode Cepheids, based on the data released from OGLE-III. We have 
applied an extinction map to correct for the extinction of these Cepheids. In addition to the VIW 
band P-L relations, we also include JHK and four Spitzer IRAC band P-L relations, derived by 
matching the OGLE-III Cepheids to the 2MASS and SAGE datasets, respectively. We also test the 
non-linearity of the Cepheid P-L relations based on extinction-corrected data. Our results (again) 
show that the LMC P-L relations are non-linear in VIJH bands and linear in KW and the four 
/i?^C bands, respectively. 
Subject headings: Cepheids — distance scale 



1. INTRODUCTION 

Recently, the Optical Gravitational Lensing Exper- 
iment (OGLE) team released a catalog of the Large 
Magellanic Cloud (LMC) Cepheids from its third phase 

• of observation ( herea fter OGLE-III), as described in 

I iSoszvhski et al.l (|2008j ). In their paper, the Cepheid 
period-luminosity (P-L) relations in the V and / 
bands were derived from this OGLE-III catalog with- 

' out any extinction correction, and it was left for other 
researchers to correct for the extinction. On the 

■ other hand, the current optical band LMC P-L re- 
lations available in the literature (see, for example, 
Udalski et al.lll999at ISandage et"an 120041: iKanbur et all 
200l IKanbur fc Ngeowl l2006t iFouaue et al.l 12001 are 
mostly derived from the OGLE-III p redecessor, the 
OGLE -II catalog (U dalski eTap [19993) ■ In the JHK 
band, iPersson et al.i ((2004) derived the P-L relation 
based on the complete ligh t cur ves of 92 LMC Ceph eids, 
while iGroenewegenl (|200C)D and iFouaue et al] (|2007f) de- 
rived the JHK P-L relations by matching the 2MASS 
(Two-Micron All Sky Survey) point sources to the 
OGLE-II Cepheids. At longer wavelengths, the Spitzer 
IRAC band P-L relations were derived by matching 
the OGLE-II Cephe ids and t he Cepheids sample froin 
IPersson et all (|2004|| . done in iNge ow fc Kanbuil (120081 ) 
and iFreedman et all (|2008D respectively, to the SAGE 
(Surveying the Agents of a Galaxy's Evolution) archival 
data. 

In this Paper, we attempt to derive the extinction cor- 



rected VI band P-L relations by applying an extinction 
map to the OGLE-II I Cepheids (which is not done in 
ISoszvhski etalll2008f ). We also derive the JHK and 
the IRA C band P-L relations by matching the OGLE- 
III Cepheids to the 2MASS and SAGE database, re- 
spectively. In addition, one unexpected result from 
the OGLE-II based P-L relation is that the optical 
LMC P-L relation was found t o be non-linear (see, for 
example, [S andage e,iafl l2004t IKanbur fc Ngeowl 12004 
iNgeow etal . 2008, and reference therein). We will also 
test the non-linearity of the extinction corrected OGLE- 
III P-L relations using rigorous statistical tests. 

2. THE DATA 

The data for OGLE-III fundamental mode (FU) 
Cepheids were kindly provided for us by I. Soszyhski: 
this included the celestial coordinates, the mean VI 
band magnitudes, and the periods for 1848 Cepheids. 
Further details abou t this dataset can be found in 
ISoszvhski et al.l ()2008D . Throughout this Paper, we have 
assumed the classification of these FU Cepheids, as given 
in Soszyhski et al. (2008), is robust. Five Cepheids with- 
out V and / band mean magnitudes were excluded from 
the sample. We also excluded the 11 Cepheids with 
logP > 1.5 (where P is the pulsation period) that only 
have V band data. The observed magnitudes of these 
Cepheids are close to the s aturation limit of th e OGLE- 
III survey (at / '--^ 13 mag.. ISoszvhski et al]|2008l ). hence 
their photometry may be affected and/or biased. 

We also matched the OGLE-III FU Cepheids with 
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the 2MASS po int source catalog (jCutri et al.l 120031 : 
iSknitskie et"al][200l . using a search radius of 1" , in or- 
der to derive the corresponding P-L relations in the JHK 
bands. The match did not find the 2MASS counterparts 
for 71 out of the 1832 Cepheids in our sample within the 
search radius. The remaining 1761 Cepheids only have 
one matched 2MASS point source with a mean separa- 
tion of 0.108" (RMS = 0.110"). W e appl ied the pre- 
scription outlined in iSoszvnski et al.l (|2005f ) to convert 
the single-epoch 2MASS photometry to the mean JHK 
magnitudes, using the epoch of the 2MASS data, the pe- 
riod and the epoch at maximum light of the Cepheids, 
and the scaling between the / band and JHK band am- 
plitudes. 

Finally, we matched the OGL E-III FU Cepheids with 
the SAGE data (.Meixner et all .2006. ) from Spitzer ob- 
servation. The IRAC band data was taken from the 
SAGE Winter '08 I RAC Archived In contr ast to the 
SAGE data used in iNgeow k. Kanbud ()2008l ). this new 
SAGE data contains the photometry from both Epoch 

1 and 2 observations. Using the same 1" search radius, 
the query to the SAGE Archive returned 1781 and 1774 
unique sources for the Epoch 1 and 2, respectively, and 
1759 of them have both the Epoch 1 and 2 photome- 
try. The mean separations for the matched Epoch 1 and 

2 sources are 0.222" (RMS = 0.114") and 0.163" (RMS 
= 0.147"), respectively. There are a total of 1796 sources 
after combining the data from the two epochs. We took 
the averaged values from the two epochs for the matched 
sources, or kept the single epoch magnitudes if either one 
of the epochs is available. 

3. THE EXTINCTION CORRECTED P-L RELATIONS 

There are various ways to estimate the extinction 
for Cepheids (for e xample, see the introduction in 
iKovtvukh et al.ll2008f). In this Paper, w e adopted the ex- 
tinction map from lZaritskv et aTl ()2004L hereafter Z04) as 
a mean to estimate the extinction for individual OGLE- 
III FU Cepheids, because the Z04 extinction map cov- 
ers almost the entir e LMC. The extinction map from 
llmara fc Blit j ()2007[ ) covers a similar region as in Z04 
extinction map, but their extinction map is not pub- 
licly available. The OGLE team has also published the 
LMC extinction ma p based on the OGLE-II red-clump 
stars ( Udalski ct al| ll999bD . which was later refined by 
ISubramaniam (|200a ). however the OGLE-II extinction 
map only covers the central bar region in the LMC. 
Comparison of the different extinction maps and var- 
ious methods is beyond the scope of this Paper, and 
will b e presented in future papers. As in I Ngeow et al.l 
(|2005h . we selected the extinction values derived from 
the cool stars when applying the Z04 extinction map. 
For a given Cepheid location, the map returns the V 
band extinction Ay- For the three Cepheids that re- 
turned a null value of Ay, we used the mean extinction 
value of E{B -V) = 0.10 appropriate for the LMC. The 
averaged Ay for the OGLE-III Cepheids is 0.429, which 
translates to a mean E{B — V) of 0.129. For extinction 
in other bands, we scaled the Ay values with R\/Ry, 
where R is the total-to-selective extinction coefficient. 
The adopted values for R are -Ry,7,,/,_f/,A',3. 6,4. 5, 5.8,8.0 = 

^ Using the IRSA's Gator Catalog Q uery at 
|http : //irsa. Ipac . caltech. edu/applicat ions/Gat or/ | 
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Fig. 1. — The V , I and W band P-L relations derived from the 
OGLE-III fundamental mode Cepheids. The V and I band mean 
magnitudes have been corrected for extinction using the extinction 
map from Z04. The filled circles are the Cepheids used to derive 
the P-L relations, and the open squares are the rejected outliers. 
The crosses are the excluded 11 long period Cepheids as described 
in the text. 
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Fig. 2. — The J, H and K band P-L relations derived from 
matching the OGLE-III fundamental mode Cepheids with 2MASS 
point source catalog. The JHK band magnitudes have been cor- 
rected for random phase and extinction as described in the text. 
The symbols are same as in Figure [T] 

{3.24, 1.96, 0.95, 0.59, 0.39, 0.17, 0.12 0.08, 0.05) 
The va lues in VI band are taken from lUdalski et al.l 
()1999bD . while for othe r band s, the values are derived 
from the ICardelli et al.l (|1989( ) extinction law. In addi- 
tion, we also derived the reddening-free Wesenheit func- 
ti on W = I —1.5 5 (V — I ), with the same expression given 
in lUdalski et al.l ()1999af ). in this Paper. 

Outliers of the P-L relations were removed using 
an iterative 2.5g cl ipping algorithm (for example, see 
lUdalski et "al]ll999af l. where a is the RMS derived from 
fitting the P-L relation to the data in a given band. The 
resulting P-L relations are presented in Figure [Tl [2] and 
[3l respectively. For the 5.8/im and 8.0/ini P-L relations. 
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TABLE 2 

Comparison of Various P-L Relations. 



log(P) 

Fig. 3. — The I R AC hand P-L relations derived from matching 
the OGLE-III fundamental mode Cepheids with SAGE Archival 
data. The symbols are same as in Figure [T] The dashed and 
dotted lines on the lower two panels are the detection limit and 
the adopted period cut in those bands, respectively. 



TABLE 1 

The LMC P-L Relation in Various Bands. 



Band 



Slope 



Zero-Point 



V 


-2.769 ±0.023 


17.115 ±0.015 


0.204 


I 


-2.961 ±0.015 


16.629 ±0.010 


0.131 


J 


-3.115 ±0.014 


16.293 ±0.009 


0.121 


H 


-3.206 ±0.013 


16.063 ±0.008 


0.107 


K 


-3.194 ±0.015 


15.996 ±0.010 


0.124 


3.6/xm 


-3.253 ±0.010 


15.967 ±0.006 


0.087 


4.5/im 


-3.214 ±0.010 


15.930 ±0.006 


0.089 


5.8/im 


-3.182 ±0.020 


15.873 ±0.015 


0.122 


8.0/^m 


-3.197 ±0.036 


15.879 ± 0.034 


0.139 


W 


-3.313 ±0.008 


15.892 ±0.005 


0.069 



we an applied additional period cut at the short period 
end. This is because the magnitudes for the short pe- 
riod Cepheids are approaching the detection limits at 
~ 15. 2maK. and ~ 14.2m ag. in the 5.8/j,m and S.O^m 
band (jMeixner et al.l[206 d). respectively. Further, the 
8.0/im P-L relation displays a flattening for Cepheids 
with logP < 0.7, which is probably due to the larger 
photometric errors when approaching the limiting mag- 
nitudes. A period cut is needed to remove the bias due 
to this effect. We adopt a period cut at logPcut 0.536 
and ^ 0.702 for the 5.8/im and 8.0/^m P-L relations, re- 
spectively. These period cuts were estimated from a non- 
linear parametric estimation procedure (same procedure 
to estimate the break period in the next section). The 
slopes become shallower without these period cuts. The 
linear version of the P-L relations in all bands is sum- 
marized in Table [TJ Note that we obtai ned an (almost) 
identical W band P-L relation as given in lSoszvhski et all 
(I2OOI . 

3.1. Comparison with the Published Results 
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Zero-Point 


N 


|T| 


t 






V Band 










_2 7f;q 4- n no*? 


1 7 1 1 K -1- f) ni t; 

-L 1 . J. iij _1_ \J ,\J ±ij 


1675 






2 


— 9 7fin -1- 0*^1 


1 7 049 -1- 091 


649 


208 


836 


Q 




9 77Q _|_ n 0*^1 


1 7 Ofifi -1- n 091 




232 


0^17 
U.Ol ( 


4 


—9 7*^4 -1- n n9Q 


1 7 n'^9 -1- 007 


716 


0.938 


0.348 














\ 




1 fi9Q -1- 01 
±u.u.^y ^ u.u±u 


1640 






2 


— 9 QR9 -1- 091 




658 


0.036 


0.972 


3 


—9 Q7Q -1- 091 


1 i^QA -1- 01 4 


662 


0.636 


0.525 


4 


_o QK7 _|_ n non 




692 


160 


0.873 






W Rand 








\ 


_Q Q1 Q _|_ n 00)^ 
— O.OJ-O ^ u.uuo 


1 ^ SQ9 -1- OO'^ 


1501 






2 


Q 977 -U 01 4 




690 


2 321 


020 


3 


—3.309 di 0.011 


15.875 zb 0.007 


671 


0.269 


0.788 


4 


— "^ '?90 -1- 01 1 


1 ^ S)^n -1- 00"^ 


686 


0.517 


0.605 






U JJOjliU. 








\ 


_Q llf;-|-n014 


1 90*^ -1- OOQ 
±u..^yo ^ u.uuy 


1586 






4 


_'3 1 QQ _|_ n r\o(< 
Kj-i-KiZj _i_ Kj.yj^yj 


1 fi 97'^ + D OOfi 

-LU . 1 vJ _1_ U.UUU 


529 


0.836 


0.403 














\ 


— 90^ -1- 01 


1 (\ nfi*^ -1- n 00^^ 

±u.uuo ^ u.uuo 


1561 






4 


_'3 oQy -1- 094 


1 n'^9 -1- OO'^ 

XU.UiJ^ ^ U.UUtJ 


529 


1.186 


0.236 






_f\. JJdliU. 








\ 


_Q 1 Q4 + 01 'S 


1 K qof! -1- f) ni n 


1554 






4 


— 99)^ -1- 09^^ 


1 ^ Q^^Q -1- n nofi 

-Ltj.yoy ^ u.uuu 


529 


1 120 


263 






3.6/jm Band 
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-3.253 ±0.010 


15.967 ± 0.006 


1617 






5 


-3.263 ±0.016 


15.945 ±0.012 


628 


0.566 


0.572 






4.5/im Band 








1 


-3.214 ±0.010 


15.930 ± 0.006 


1633 






5 


-3.221 ±0.017 


15.927 ±0.012 


635 


0.377 


0.707 






5.8/^m Band 








1 


-3.182 ±0.020 


15.873 ±0.015 


931 






5 


-3.173 ±0.028 


15.850 ±0.022 


561 


0.272 


0.785 






8.0/im Band 








1 


-3.197 ±0.036 


15.879 ±0.034 


401 






5 


-3.091 ±0.039 


15.684 ±0.036 


319 


1.954 


0.051 



Note . 
Ill999al): 



Reference: (1) this work; (2) 



(3) same as (2) but with the updated version given in 

|f tp: //sirius . astrouw . edu . p i/ ogle/ogle2/var_stars/lmc/cep/catalog/] 
README. PL; (4) lFouaue et al] | [20071) : (5) INgeow &: Kanburl ll2008f) . 



In Table [H we compare our multi-band P-L rela- 
tions from Table [1] with the published optical VIW 
band P- L relations de r ived from the OGLE -II Cepheids 
fUdalski et all Il999at iFouaue et aU l2007l). as well as 
the JHK band P-L relations (iFouaue et all l2007f) and 
the IRAC band P-L relations (|Ngeow fc KanbuH I2OO8I) 
that matched the OGLE-II Cepheids to the 2MASS and 
SAGE (Epoch 1) data, respectively. In this comparison, 
our philosophy is that it is not sufhcient to simply com- 
pare whether the P-L slope and/or zero-point obtained 
by one sample is within one or two standard deviations 
of the P-L slope and/or zero-point from another sample: 
this is a necessary but not sufficient condition. What 
is needed is a simultaneous comparison. An analytical 
argument showing this has been stated in iNgeow et al.l 
(|2008?). 

We applied the standard t-test (for example, see 
iBethea et al.l [19951 ) to test the published P-L slopes in 
Table [21 under the null hypothesis that the slope derived 
in this Paper, for a given band, is consistent with the 
published slopes. The T- values calculated from the t-test 
incorporate the standard deviations of the slopes and the 
sample variance in the two samples under testing. The 
corresponding i(Q/2.iy)"Values were then evaluated from 
the t distribution with v — Ni + N2 — ^ degree of free- 
dom, where A^i and N2 are the number of Cepheids in 
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the two samples. In our test, we adopted a constant 
significance level of a = 0.05. The values of T and t for 
each pair of the slopes were given in the last two columns 
of Table [2l The null hypothesis of equal slopes can be 
rejected if \T\ > t. 

In general, both agreements and disagreements were 
found when comparing the slopes to the published re- 
sults. The disagreements include the VJH K band slopes 
when compared to those published in iFouque et al.l 
(|2007f ). t he / band slope when compared to the updated 
version of lUdalski efall (|1999aD . the W band slope whe n 
compared to the slope given in lUdalski et al"] (Il999af ). 
and the S.O^tm band slope w hen compared it to the slope 
found in lNgeow fc Kanburl {2008). However, agreements 
were also found for other VIW slopes in the same band. 
In all cases, the slopes that are in disagreement are either 
shallower (in VW band) or steeper (in / band) than the 
rest of the slopes in the same bands. In the JHK band, 
the slopes fro m Table [1] were shallow er than the counter- 
parts given in iF ouque et al.l (|2007f l. Since both sets of 
P-L relations matched the LMC Cepheids to the 2MASS 
data and applied the same transformation of the single 
epoch photometry to the mean magnitudes (as given in 
ISoszvhski et al.|[2005f ). the discrepancy of the slopes is 
due to the different sample size used in both studies 
(there are about t hree times more Cep heids in our sam- 
ples than those in IFouque et al.l 120071 )^. The slopes in 
3.6/j,m, 4.5Atm and 5.8jum P-L rel ations agree with those 
given m UN geow k Kanbud (|2008f l. though the result for 
3.6/im band is marginal. The discrepancy for the 8.0/im 
band P-L relation is mainly due to the additional period 
cut applied in this Paper. 

The difference of the zero-points in the VI band P- 
L relations between Table [T] and the published results 
is mainly due to the different extinction maps used in 
der iving the P-L rela tions. The P-L rel ations given in 
lUdalski et all dUlli) and IFouque et all (|2007l ) were de- 
rived using the OG LE-II extinction m ap, with a mean 
E{B -V) ^ 0.147 mdalski et an fT999a). This value is 
higher than the mean E{B — V)= 0.129 found in this 
paper from the Z0 4: extinction map (see a similar find- 
ing , for example, in ISubramaniamll2005l : iFingerhut et al.l 
l2007l ). Hence the derived P-L zero-points will be fainter 
by ~ 0.06mag. and ^ O.OSmag. in the V and I band, re- 
spectively, to those OGLE-II based P-L relations. This 
difference is indeed seen i n Table [H (except fo r the / 
band P-L zero- point from lUdalski et cd] Il999a[ ). The 
difference of the P-L zero-points in other bands, that 
less affected by extinction, may be due to the different 
sample size used in d eriving the P-L relations (see also 
ISoszvnski etalll2008f ). 

4. TESTING FOR NON-LINEAR P-L RELATIONS 

As in previous work, the non-linearity of the P-L re- 
lation is referred to the two P - L relations separated at 
10 davs (ISandaee et al.| 120041 : iKanbur fc Ngeowl l2004t 
iNgeow et al]l2005ll2008( ). We have applied a non-linear 
parametric estimation procedure available in the SAS 

^ If restricted to ~ 610 OGLE-III Cepheids that are com- 
mon to OGLE-II, the P-L slopes we obtained in JHK band are 
-3.135 ± 0.022, -3.228 ± 0.020 and -3. 213 ± 0.024, r espectively. 
When compared to the P-L relations from lFouque et all |2007D, the 
corresponding T- and t- values are \Tj h k\ = {0.118, 0.290, 0.409} 
and = {0.906,0.772,0.683}. 



package, which utilizes the Levenberg - Marquardt al- 
gorithm, to estimate the break period from the / band 
data (because the / band light curves are well sampled 
from large number of observations) . The estimated break 
period is located at logPbrcak = 1.006 (with the 95% con- 
fidence limits given by ±0.170), which is fully consistent 
to the adopted 10 day period in defining the non-linearity 
of the P-L relation. 

The fitted long (logP > 1.0) and short period P-L 
relations are summarized in Table [31 At first glance, 
the difference of the slopes between the long and short 
period P-L relations in V, K, 3.6/Ltm, 4.5/Ltm, 5.8/im and 
W band seems to be within ~ 1.5cr. This may suggest 
these P-L relations are linear. In contrast, the I, J, H 
and 8.0/xm band P-L relations show a difference of the 
slopes at ~ 2. Oct level or larger. We emphasize again 
that testing for non-linearity by simply seeing if the long 
and short period slopes are within a certain number of 
standard deviations of each other is a necessary but not 
a sufficient condition: a simultaneous comparison using 
statistical test is needed. Thus a P-L relation can be 
non-linear even if the short and long period slopes are 
wi thin 1 to ~ 2ct of e ach other as an alytically outlined 
inlNgeow et al.l (|2008f) . Furthermore. iNgeow fc Kanbuil 
()2006[) also argued that statistical tests are needed to 
detect the non-linear P-L relation. 

In this Paper, we only use the F-test, though sev- 
eral other statistical tests have b een applied to the 
OGLE-II based P - L relations before (iKanbur et"all l2007l: 
iKoen et al.l 120071 : INgeow et ar '2008! ). De tails for the 
F-test can be fou nd in Kanbur fc Ngeowl ([2004) and 
INgeow et al.l ()2005D . and will not be repeated here. For 
a large number of data points, F ~ 3.0 at 95% confi- 
dence level. Hence, the underlying P-L relation is sta- 
tistically consistent with non-linearity if F > 3 (that is, 
the null hypothesis of linear regression can be rejected). 
The choice of a 95% confidence level is fairly standard 
and consistent with our previous work. 

The robustness of the F-test has a large literature. 
Its assumptions are homoskedasticity, normality of er- 
rors and independent observations. The last assumption 
regarding independent observations means that observa- 
tions of one Cepheid have no b earing on subsequ e nt ob - 
servations of another Cepheid. IKanbur fc Ngeowl (|2004l ) 
checked assumptions of normality and homoskedasticity 
using OGLE-II Cepheids data and found them to be rea- 
sonable assumptions. Homoskedasticity essentially re- 
quires there to be not too many outlier s in the data, 
which have been removed in our samples. lAli fc Sharmal 
()1996l ) suggest the F-test is robust to small departures 
from normality and our large sample size, both with 
the OGLE-II and OGLE-III Cepheids, further adds to 
this robustness. For example, INgeow fc Kanbuj ()2006D 
appended data from different sources to increase the 
OGLE-II sample size and found the F-test produced con- 
sistent results. 

The F-test results were given in the last two columns 
of Table |3l From this table. The F-test results clearly 
suggest that the IJH band P-L relations are non-linear, 
in contrast the K band P-L relation and the 3.6^m to 
5.8/im P-L relations are linear. These results support 
the previous finding of the non-linear/linear LMC P- 
L relations from the OGLE-II and MACHO Cepheids 
(jNgeow et al.ll2005HKa"nbur fc Ngeowll2006l : iNgeow et all 
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TABLE 3 

F-Test Results of the P-L Relations. 



Band 


Slope5 


Zero-Point5 




Ns 


Slope 


Zero-Point/^ 




Nl 


F 


p{F) 


V 


-2.823 ±0.031 


17.143 ±0.018 


0.202 


1566 


-2.746 ±0.165 


17.122 ±0.195 


0.230 


109 


3.18 


0.042 


I 


-3.004 ±0.020 


16.651 ±0.012 


0.130 


1553 


-2.775 ±0.111 


16.440 ±0.129 


0.132 


96 


5.77 


0.003 


J 


-3.150 ±0.019 


16.312 ±0.011 


0.119 


1486 


-2.909 ±0.120 


16.075 ±0.139 


0.142 


100 


5.24 


0.005 


H 


-3.246 ±0.017 


16.085 ±0.010 


0.106 


1461 


-2.989 ±0.096 


15.832 ±0.113 


0.117 


100 


7.90 


0.000 


K 


-3.212 ±0.021 


16.006 ±0.013 


0.124 


1449 


-3.057 ±0.101 


15.845 ±0.117 


0.125 


105 


1.57 


0.209 


3.6/im 


-3.254 ±0.013 


15.967 ±0.008 


0.086 


1509 


-3.233 ±0.068 


15.944 ±0.082 


0.100 


108 


0.06 


0.946 


4.5/im 


-3.228 ±0.014 


15.938 ±0.008 


0.088 


1524 


-3.185 ±0.071 


15.903 ±0.085 


0.104 


109 


1.05 


0.352 


5.8/im 


-3.164 ±0.035 


15.861 ±0.025 


0.122 


813 


-3.275 ±0.083 


15.984 ± 0.099 


0.125 


118 


0.78 


0.460 


8.0/im 


-3.036 ±0.100 


15.746 ± 0.083 


0.142 


286 


-3.348 ±0.086 


16.058 ±0.102 


0.127 


115 


2.70 


0.069 


W 


-3.329 ±0.010 


15.900 ±0.006 


0.068 


1585 


-3.338 ±0.084 


15.934 ± 0.096 


0.084 


84 


2.90 


0.055 



Note. — cr is the dispersion of the P-L relation. Subscripts s s.nd l refer to the short (log P < 1.0) and long period Cepheids, 
respectively. The p{F) is the probability under the null hypothesis of single line regression, for the given F value and the degree 
of freedom. 



120081 iNgeow fc KanbuH l2008f) . and extend the non- 
linear/linear results further to the Cepheids from the 
OGLE-III sample. The F-test resuhs on the V, W and 
8.0/im P-L relations are worth a further discussion. 

In our previous work, the LMC V band P-L relation 
shows a hig hly significa nt f -test result for non-linearity 
(INgeow et a l. 2005: Kanbur fc Ngeowll2006t iNgeow et ail 
|2008[) . with a large F- value and a p- value that is signif- 
icantly smaller than 0.05. The F-value listed in Table [3] 
is 3.18, which appears to be close to the F — 3.0 thresh- 
old at 95% confidence leveP However, hypothesis testing 
works with p-values for the chosen statistics under the 
null hypothesis. Hence, the probability of having an F- 
value as large as 3.18 is smaller than 0.05. Therefore, the 
V band P-L relation does show evidence of non-linearity 
from the F-test. 

For the Wesenheit function, INgeow fc Kanbuil ()2005l ) 
and INgeow et all (|2008f) found that the LMC Wesenheit 
function, using the OGLE-II data, derived from the op- 
tical band is linear. The F = 2.90 (with p — 0.055) from 
Table |3| suggests the W band P-L relation found in this 
paper is also linear. However, if the p = 0.042 for the 
computed V band F statistics in Table [3] is considered 
to be non-linear, then the same attention must be given 
to the fact that the p = 0.055 for W band P-L relation 
which is close to being considered linear. Our explana- 
tion for this has been that the non-linearities in the P-L 
and period-color (P-C) relation s cancel out in produc- 
ing a linear Wese n heit function fNgeow fc KanbuJ l2005l : 
iKoen et al]l2007| ). iKoen et all (|2007l ) also suggested that 
with a larger sample and more accurate data, the "can- 
cellation" of non-linearities in the PL/PC relation will 
not necessarily be exact leading to a non-linear Wesen- 
h eit function. 

INgeow fc Kanbuil (|2008f) found the 8.0/im P-L relation 
is non-linear, in contrast to the linear F-test result given 
in Table |3l This discrepancy is due to the period cut 
applied at a rath er large period (note th at no period 
cut was appfied in INgeow fc Kanbuil 120081 ). In fact, the 
8.0/zm P-L relation becomes non-linear if logPcut < 0.69 
is adopte d. The n on-linear 8.0/im P-L relation found in 
[Ngcow & KanbuH (2008) could be caused by the incom- 
pleteness of the data at the faint end when approaching 
the detection limit. Another possibility is that mass- loss 

^ If the confidence level is lowered, say to 90%, then the cor- 
responding _F-value will be decreased to ~ 2.3. In this case the 
F-test result for the V band P-L relation clearly shows that it is 
non-linear. 



may plays some role in longer wavele ngth P-L relat ions to 
make the P-L relations non- linear (jNeilson et al.l[2008. ). 
Detailed study of the effect of mass-loss on the IRAC 
band P-L relations, based on the OGLE-III Cepheids, 
is currently underway, and will be presented in a future 
paper. 

5. INCLUDING THE LONGER PERIOD CEPHEIDS 

In previous sections, we excluded the Cepheids with 
logP > 1.5 from the OGLE-III sample, which only have 
the V band data. Including these Cepheids the linear 
version of the V band P-L relation is almost identical 
to the relation given in Table [1] However, the slope of 
the long period P-L relation becomes steeper (—2.823) 
as compared to the slope given in Table |3| (—2.746) if 
these Cepheids wer e included. Previous studies (see, for 
example, [S andagc et all l2004t iKanbur fc Ngeowl 120061 : 
INgeow et al.l 120 08) found that the long period V band 
P-L relation is shallower than the short period counter- 
part. Hence the steeper slope may suggest these longer 
period V band Cepheid data should not be included in 
the sample. 

It is possible to increase the number of long pe- 
riod Cepheids and extend the period coverage beyond 
logP = 1.5, by including the additional Cepheids from 
the literature. This has been done, for example, in 
ISandage et al.l ()2004l ). However, one criticism has been 
that the published long period Cepheid data may not 
be i n the same photome tric system as in the OGLE sur- 
vey (jFouque et al.ll2007l ). The on-going "OGLE Shallow 
Survey of the LMC Cepheids" is aimed to observe the 
longer period LMC Cepheids in the same photometry 
system as in OGLE-III. Hence, it is better to wait for 
the data released from this survey before appending ad- 
ditional longer period Cepheids to our sample. 

6. CONCLUSION 

In this Paper we derive the VIJHKW and the four 
ffl^Cband LMC P-L relations based on the OGLE-III 
Cepheids. The Cepheid data were corrected for extinc- 
tion by applying the Z04 extinction map. We matched 
the OGLE-III Cepheids to the 2MASS point sources, and 
the single-epoch 2MASS photometry was transformed 
to the mean magnitudes. For the IRAC band data, a 
straight average from the SAGE Epoch 1 and 2 archival 
data were adopted, which s hould closer to the mean mag- 
nitudes than those used in INgeow fc Kanb ur (200^). 

We compare our P-L relations to the published results 
using the i-test. For the published optical P-L relations 
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that used in the comparison, 2/3 of them show a good 
agreement on the P-L slopes to our results. For the JHK 
and IRAC band, the i-test results show that our P-L 
slopes disagree and agree (except for the 8.0/im band 
P-L slope) to the published results, respectively. 

We further test the non-linearity of the P-L relations 
using the i^-test to the OGLE-III Cepheids, and found 
that the LMC VIJH band P-L relations are non-linear 
but linear in KW and the four /i?^C bands, respectively. 
It has to be cautious that the 5.8/im and 8.0/xm P-L rela- 
tions may be affected by incomplete bias at the faint end 
and/or mass-loss. It is worth pointing out that both 
MACHO and OGLE-III samples show the same non- 
linear/linear results in V JHK bands, where both sam- 
ples have relatively large number of Cepheids and cover 
the similar region in the LMC. 
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